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We investigate the statistics of fluctuations in a classical stochastic network of nodes joined by 
connectors. The nodes carry generalized charge that may be randomly transferred from one node 
to another. Our goal is to find the time evolution of the probability distribution of charges in the 
network. The building blocks of our theoretical approach are (1) known probability distributions for 
the connector currents, (2) physical constraints such as local charge conservation, and (3) a time- 
scale separation between the slow charge dynamics of the nodes and the fast current fluctuations 
of the connectors. We integrate out fast current fluctuations and derive a stochastic path integral 
representation of the evolution operator for the slow charges. The statistics of charge fluctuations 
may be found from the saddle-point approximation of the action. Once the probability distributions 
on the discrete network have been studied, the continuum limit is taken to obtain a statistical field 
theory. We find a correspondence between the diffusive field theory and a Langevin equation with 
Gaussian noise sources, leading nevertheless to non-trivial fluctuation statistics. To complete our 
theory, we demonstrate that the cascade diagrammatics, recently introduced by Nagaev, naturally 
follows from the stochastic path integral. By generalizing the principle of minimal correlations, 
we extend the diagrammatics to calculate current correlation functions for an arbitrary network. 
One primary application of this formalism is that of full counting statistics (FCS), the motivation 
for why it was developed in the first place. We stress however, that the formalism is suitable for 
general classical stochastic problems as an alternative approach to the traditional master equation 
or Doi-Peliti technique. The formalism is illustrated with several examples: both instantaneous and 
time averaged charge fluctuation statistics in a mesoscopic chaotic cavity, as well as the FCS and 
new results for a generalized diffusive wire. 



PACS numbers: 73.23.b, 02.50.r, 05.40.a, 72.70.+m 
I. INTRODUCTION 

Consider an exclusive night-club with a long line at the 
entrance. A bouncer is at the front of the line to keep 
out the rif-raf. At every time step, a person is accepted 
inside the club with probability p, or rejected with prob- 
ability 1 — p. Inside the club, people stay for a while and 
eventually leave. At every time step, the probability a 
person leaves is q. We want to answer a question such 
as "what is the probability that Q people leave the club 
after t time steps?" . 

Assuming that p and q remain constant, the situation 
is simple and we can easily solve the relevant probabilis- 
tic problem. However, in realistic situations this rarely 
happens: the management wants to make money. If the 
club is almost empty, they instruct the bouncer to be 
less discriminating, while if the club is almost full, the 
bouncer is to be more discriminating. Thus, p becomes a 
function of the number of people in the club. People will 
be more likely to leave if the club is very crowded, so q is 
also a function of the number of people inside the club. 
The problem posed now is much more difficult because 
of the presence of feed-back: the elementary processes 
change in response to the cumulative effect of what they 
have accomplished in the past. 

This simple example captures all the basic features of 
the problems we wish to consider. Although the exam- 
ple was given with people, the actors in the probability 
game may be any quantity such as charge, energy, heat 
or particles, which we will refer to simply as generalized 



charge. Similarly, the night club can be a mesoscopic 
chaotic cavityji a birth-death process^ a biological mem- 
brane channel^ etc. 

Historically, general stochastic problems are solved 
with the master equation. The time rate of change of 
the probability to be in a particular state is given in 
terms of transition rates to other states. This approach 
has had great success and leads naturally to the Fokker- 
Planck and Langevin equations*^ However, once the mas- 
ter equation is given, the solution is often quite difficult 
to obtain. 

This paper takes a different approach. Rather than 
beginning with a master equation describing the prob- 
ability of all processes happening in a unit of time, we 
make several assumptions from which we can reformulate 
the problem. Although these assumptions limit the ap- 
plicability of the theory, when they apply, the problems 
are much easier to solve. The assumptions are: 

• The system we are interested in is a composite sys- 
tem made out of constituent parts. In the night 
club example, the system is made up of three phys- 
ical regions: outside the front door, the interior of 
the club, and outside the back door. The decom- 
position of a larger system into smaller interacting 
parts is only meaningful for us if there is a separa- 
tion of time scales. This means that the charge in- 
side the constituent parts changes on a slower time 
scale than the fluctuations at the boundaries. In 
the night club example, this simply means that the 
average time a person spends in the club will be 
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much longer than the typical time needed to enter 
the door. 

• Taken alone, the parts of the composite system 
have a finite number of simple properties or pa- 
rameters. The only property of the night club that 
was relevant for the problem was the total number 
of people in it at any given time. The important 
element of the line out in front is that it never runs 
out. All other details are irrelevant. 

• In the limit where all parts of the network are very 
large (so that the elementary transport processes do 
not affect themselves in the short- run), the trans- 
port probability distributions between elements are 
known. In the night club example, the probability 
of getting Q people through the front door after t 
time steps (given a constant, large number of peo- 
ple inside) is easy to find, because we have assumed 
that the elementary probability p does not change 
from trial to trial. The transport probability distri- 
bution is simply the binomial distribution- where 
the probability p is a function of the (approximately 
unchanging) number of people inside. The back 
door distribution is obtained in the same way. 

• There are conservation laws that govern the proba- 
bilistic processes. No matter what probability dis- 
tributions we have, there are certain rules that must 
be obeyed. The net number of people that enter, 
stay, and leave the club must be a constant. This 
means that the time rate of change of the club's oc- 
cupancy is given by the people-current in minus the 
people-current out. The people in the line outside 
are a special case. There is in principle always a 
replacement, so moving one person inside the club 
doesn't affect the properties of the line. 

Now, the strategy is to use this information as the 
starting point to find transport statistics for the com- 
bined interacting system. The main result derived is 
a path integral expression for the conditional probabil- 
ity (taking conservation laws into account) for starting 
and ending with a given amount of charge at each loca- 
tion after some time has passed. From this conditional 
probability, specific quantities such as transport statis- 
tics through the system, fluctuation statistics of charge 
at a particular location and the like may be found. 

One primary application of this formalism is that of 
full counting statistics (FCS)^*& the motivation for why 
it was developed in the first placed FCS describes the 
fluctuations of currents in electrical conductors. It gives 
the distribution of the probability that a certain num- 
ber of electrons pass a conductor in certain amount of 
time. Mean current flow and and shot noisei corre- 
spond to the first and second cumulant of this distri- 
bution. The full distribution (defined by all cumulants) 
provides a full characterization of the transport proper- 
ties of a electrical conductor in the long time limit. In the 



past, FCS was mainly addressed with quantum mechan- 
ical tools such as the scattering theory^*^ of coherent 
conductors, the circuit theory based on Keldysh Green 
functionsjiSiiLiSii^ or the nonlinear a model^ However, 
a number of works realized that for semi-classical systems 
with a large number of conductance channels, shot noise 
may be calculated without accounting for the phase co- 
herence of the electroniA£ii£*iLi& These works treat the 
basic sources of noise quantum mechanically, but calcu- 
late the spread of the noise throughout the conductor 
classically. For specific conductors like diffusive wires 
and chaotic cavities, this idea has been extended to the 
calculation of third and fourth cumulants via the cas- 
cade principl e*^ ' and to the full generating function 
of FCSiZ'2i*22i2£ In the present work, we consider an 
abstract model instead of any particular example and 
develop the mathematical foundations of the proposed 
semi-classical procedure to obtain FCS. We introduce 
and investigate networks of elements with known trans- 
port statistics and show how the FCS of the entire net- 
work can be constructed systematically. 

The formalism we present is related to a different ap- 
proach in non-equilibrium statistical physics called the 
Doi-Peliti technique^ The idea is that once the basic 
master equation governing the time evolution of proba- 
bility distributions is given, it may be interpreted as a 
Schrodinger equation which may be cast into a second- 
quantized language. This quantum problem is then con- 
verted into a quantum mechanical path integral (often 
obeying bosonic or fermionic statistics) from which one 
may take the continuum limit and use a field theory 
renormalization group approach with diagrammatic per- 
turbation expansion«2S> This approach is useful in many 
situations far from equilibrium and has several parallels 
to our approach. It has been pointed out that this tech- 
nique is in some sense the classical limit of the quantum 
mechanical Keldysh formalism^ the same tool used in 
the past to calculate FCS, so this gives another connec- 
tion with the subject matter we are concerned with. 

There are several advantages of our approach. First, 
we skip the master equation step. If the probability dis- 
tributions of the connector fluctuations are given, we may 
immediately construct network distributions. Second, 
from a computational view, our formulation of the prob- 
lem is much simpler than starting from first principles 
for situations where the ingredients we need are available, 
and results are much easier to obtain than beginning with 
the master equation alone. Thirdly, our formulation also 
applies to situations where temporal transition proba- 
bilities may be large. Finally, the formalism's physical 
origin is clear, so the needed mathematical objects are 
well motivated. 

The rest of the paper is organized as follows. In Sec. 
Hll we introduce and develop the general theory. After 
reviewing elements of probability theory, we derive the 
stochastic path integral for a network of nodes as well 
as explore the relationship to the master equation and 
Doi-Peliti formalism. In Sec. IIIII the continuum limit 
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FIG. 1: An arbitrary network. Each node has charge and 
counting variables {Q Q , A a }. The nodes transfer charge via 
currents I a /3 through the connectors. The absorbed counting 
fields (Aq) are constants by definition of the absorbed charges 
Q° (see text). Each node may have an arbitrary number of 
different charge species, Q a = {Qa, Q 2 a , . . . , Q J a }. 



is taken to derive a stochastic field theory and link our 
formalism with the Langevin equation point of view. In 
Sec. IIVI we develop diagrammatics rules to calculate cu- 
mulants of the current distribution as well as current cor- 
relation functions for an arbitrary network. Sec. IVl gives 
several applications of the theory to different physical sit- 
uations. We solve the field theory for the mesoscopic wire 
and demonstrate universality in multiple dimensions as 
well as present new results for the conditional occupation 
function and probability distribution. We also consider 
the problem of charge fluctuation statistics (both instan- 
taneous and time-averaged) in a mesoscopic chaotic cav- 
ity. Sec. I VII contains our conclusions. 



II. GENERAL FORMALISM 

Once we have the basic elements of our theory (the gen- 
eralized charges), we must specify some spatial structure 
that they move around on. As we noted in the introduc- 
tion, the essential structure needed to state the problem 
are simply points we refer to as nodes, joined by con- 
nectors. This defines a network (see Fig. QJ. The state 
of each node a is described by one (effectively continu- 
ous) charge Q a r^ an d Q is the charge vector describing 
the charge state of the network. The node's state may 
be changed by transport: flow of charges between nodes 
takes place via the connectors carrying currents I a p from 
node a to node 0. The variation of these charges Q a is 
given by 



Qa(t + At) — Q a (t) = ^2 Qaf: 



(1) 



where the transmitted charges Q a p(t) — J" Q At dt' I a p(t+t') 
are distributed according to P a p(Q a p(t))- The fact that 
the probabilities P a p(t) also depend on the charges Q(i) 
is one source of the difficulty of the problem. 

Assuming that the probability distributions P a p 
(which depend parametrically on the state of nodes a and 
(3) of the transmitted charges Q a p are known, we seek the 



time evolved probability distribution T(Q,i) of the set 
of charges Q for a given initial distribution T(Q,0). In 
other words, one has to find the conditional probability 
(which we refer to as the evolution operator) U (Q, Q', t) 
such that 



r(Q,t) = j dQ'[/(Q,Q',i)r(Q',0) 



(2) 



We assume that there is a separation of time scales, tq <C 
tc, between the correlation time of current fluctuations, 
To, and the slow relaxation time of charges in the nodes, 
tc- As we will show in the next section, this separation of 
time scales allows us to derive a stochastic path integral 
representation for the evolution operator, 

U(Q f ,Qi,t) = /z>Q2>Aexp{S(Q,A)}, (3a) 

5(Q,A) = / dt'[-iA ■ Q 



+ (l/2)53^(Q,A a -A /9 )], (3b) 

a/3 

where the vector A has components X a : node variables 
conjugated to the Q a that impose charge conservation in 
the network. 

In the following, we define the functions H a p as the 
generating functions of the fast currents between nodes a 
and j3. On the time scale At ^> tq, the currents through 
isolated connectors are Markovian, so that all cumu- 
lants (irreducible correlators which are denoted by dou- 
ble angle brackets) of the transmitted charge (((Q a p) n )) 
are linear in At. Following the standard notation in 
mesoscopic physics^ we define the current cumulants 
(((lap) 71 )) as the coefficients in 



(((Q a p) n )) = At(((I a p) n )), 



(4) 



where the tilde symbol has been introduced to distinguish 
the bare currents of each connector (the sources of noise) 
from the physical currents I a p flowing through that same 
connector when it is placed into the network. Then the 
generators H a p are defined via the equation 



«0M"» 



d n H a p(Q,\ a p) 



(id\ af: 



(5) 



A ofl =0 



and thus contain complete information about the statis- 
tics of the noise sources. The X a p [eventually to be re- 
placed with with A Q — Xp in Eq. |3|] is the generating 
variable for the current I a p. The notion of current cu- 
mulants is useful because they are the time independent 
objects, and thus have a time independent generators, 
Eq. (J5J. The generators H a p(Q, X a p) depend in gen- 
eral on the full vector Q and not just on the generalized 
charges of the neighboring nodes Q a and Qp. This may 
serve to incorporate long range interactions between dis- 
tant nodes. 
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The charge Q a p transferee! through the connectors 
[characterized by Eqs. 1] lloll ] may be discrete. However, 
the charge in the nodes Q a is treated as an effectively 
continuous variable in Eqs. 11 IMP . This is justified if many 
charges in the node participate in transport. Formally, 
this limit allows a saddle-point evaluation of the propa- 
gator JjJHJ. 



A. Derivation of the Path Integral. 

To derive the path integral Eq. we follow the usual 
procedure^ and first discretize time, t = nAt to derive 
an expression for U that is valid for propagation over one 
time step At. Because of the separation of time scales 
T o *C tc, we can consider At as an intermediate time 
scale, 



T < At < t c . 



(6) 



The left inequality, To <C At, implies that the transmitted 
charges Q a p are Markovian4 This means that charges 
transmitted in separate time intervals are uncorrelated 
with each other. While it is not necessary to specify the 
source of the current correlation in the general formula- 
tion, it is worth noting two examples. In a mesoscopic 
point contact, the correlation time To has the interpreta- 
tion of the time taken by an electron wavepacket to pass 
the point contact. In chemical dynamics, it could be the 
time taken for a long molecule in solution to traverse a 
filter. 

In a time At, the probability that charge Q a p is trans- 
mitted between nodes a and (3 can be written as the 
Fourier transform of the exponential of a generating func- 
tion S a p: 

exp{-iX a pQ a p + S a p{X a p)} . 
The definition of the cumulant of transmitted charge is 



«(Qa/3) n » = 



d n S a p(\ a p) 



{id\ af: 



(8) 



The Markovian assumption implies that the probability 
of transmitting charge Q a p in time At followed by charge 
Q' a p in time At' through any connector is given by the 
product of independent probability distributions. This 
implies that the probability of transmitting charge Q a p 
in time At + At' may be calculated by finding all ways 
of independently transferring charge Q' a g in the first step 
and Q a p — Q' a p in the second step, 

P(Q a p, At + At') = J dQ' Q pP(Q a p - Q' a p, At') 

xP(Q' a p,At), (9) 

which takes the form of a convolution of probabili- 
ties. Applying a Fourier transform to both sides of 



Eq. © with argument X a p decouples the convolution into 
product of the two Fourier transformed distributions. 
Eq. © implies S a p(At + At' , X a p) = S a p(At,X a p) + 
S a p(At', X a p). It then immediately follows that the gen- 
crating function must be linear in time. Therefore, a time 
independent H a p may be introduced: S a p = AtH a p. 
The linear dependence of S a p on time implies that all 
charge cumulants JSJl will be proportional to time. There- 
fore, we define the time independent current cumulants, 
Eq. ©. 

Different connectors are clearly uncorrelated for At <C 
tc , which indicates that the total probability distribution 
of transmitted charges is a product of the independent 
probabilities in each connector: 30 



P\{Q a ?}] = J] Pa(3[Q a 0,At}. 



(10) 



Thus far, the analysis is only valid for times much smaller 
than tc ■ For this case, the charges in the nodes will only 
slightly change. Since we wish to consider longer times, 
we need to take into account the fact that charge trans- 
fer between different nodes will be correlated as charge 
piles up inside the nodes. This may be accounted for 
by imposing charge conservation Eq. (y during the time 
interval with a delta function, 



S{Qc 



= f ^exp{-i\ a [Q a -Q' a -J2Q a p}}. (11) 

Here, Q' a is the charge in the node before the time interval 
while Q a is the charge accumulated in the node after the 
time interval is over. In Eq. i|ll|) . X a (referred to as a 
counting variable) plays the role of a Lagrange multiplier. 
The propagator is obtained by multiplying the constraint 
(|llfl and the independent probability distribution I|1U|) . 
Representing the probabilities in their Fourier form Q 
then yields 

a 77 a>p J n 

S = ~»y^ X a {Qa - Qq - ^ Qa/3) 
a p 

+ ^HAa /3 Qa/3 + A^a/ 3 (Q',A Q/3 )]. (12) 

The full propagator U(Q, Q', Q a p, At) still keeps track 
of each individual connector contribution Q a p- We now 
integrate out the fast fluctuations to obtain the dynamics 
of the slow variables. This may be done by using the 
identity Y, a Y.p Qap = Y, a >p(^c,Qc,p + XpQp a ) and 
Q a p — ~Qp a - The integration over Q a p gives a delta 
function of argument X a p — {X a — Xp), so that the X a p 



5 



integrals may be trivially done. We obtain 
U(Q, Q', At) = J] / ^ exp { - i£ UQ a - Q' a ) 

a a 

+At^ff Q/3 (Q',A Q -A^)}. (13) 

a>p 

This is the general result for the one step propagator. If 
any two nodes are unconnected, H a p is zero. 

An important comment is in order: because H a p 
changes slightly over the time period, which in turn af- 
fects the probability of transmitting charge through the 
contacts, it is not clear at what part of the time step 
H a p should be evaluated. This ambiguity exists because 
our theory is not microscopic. Rather, it takes the mi- 
croscopic noise generators as an input. This ambiguity 
gives the freedom of stochastic quantization^! The same 
problem also occurs in quantum mechanical path inte- 
grals, and its source there is an ambiguity in operator 
ordering^ As we are interested in the large transport- 



with 

H(Qk,Ak) = ^2 H a fi [Q a ,k',^a,k — <\g,fe] , (14b) 

a>/3 

where we have introduced the notations e?Qfc = Y[ a dQa,k 
and dA k = ]^[ Q (dA Qj fc/27r). We are now in a position to 
take the continuous time limit. Writing Q^+i Q k = 
AtQ, which is valid because the charge in any node 
changes only slightly over the time scale At, the action of 
this discrete path integral has the form S — At X)l-=i $k, 
which goes over into a time integral in the continu- 
ous limit. Using the standard path integral notation 
/PQPA = J dA rife=i / dQ k dA k , and invoking the 
symmetry H a p{X a — Xp) — Hp a (Xp — X a ) we recover 
Eq. ©. The only explicit constraint on the path inte- 
gral comes with the charge configurations at the start 
and finish, Qi and Q/. We also note that H a p depends 
on any external parameters such as voltages or chemical 
potentials driving the charge Q. 

In the simplest case of one charge and counting vari- 
able, the form of the path integral is the same as the 
(Euclidian time) path integral representation of a quan- 
tum mechanical propagator in phase space with position 
coordinate Q and momentum coordinate Xm^ The differ- 
ences with the quantum version are that the propagator 
evolves probability distributions, not amplitudes (simi- 
larly to Ref . l25f) . as well as the fact that the "Hamilto- 
nian" H = (1/2) Y. a p H ap(Q, K - Xp) is not really a 



ing charge limit, 7 > 1, and evaluate the integrals in 
leading order saddle-point approximation, this ambigu- 
ity will not affect the results^ For calculations beyond 
the large transporting charge limit, the canonical vari- 
ables Q and A need to be properly ordered, which can 
only be done with a microscopic theory. For example, 
the master equation discretized in time as discussed in 
Sec. Ill Dl requires the placement of A operators in front 
of Q operators, since the generating functions H a p of the 
transition probabilities depend on the state of the system 
at the beginning of the time period. 

To extend the propagator (|13fl to longer times t = nAt, 
we use the composition property of the evolution opera- 
tor (also known as the Chapman-Kolmogorov equation^) . 
This requires separate {Q a } integrals at each time step, 
so that for n time steps there will be n — 1 integrals over 
Q, while each of the n one-step propagators comes with 
its own A integral, A = {A Q }. Inserting our expression 
for the At step propagator Eq. (fH3)l . we find 



(14a) 



Hamiltonian, but rather a current cumulant generating 
function and therefore is not Hermitian in general. Even 
so, because of the similarity we shall refer to H as the 
Hamiltonian from now on. 



B. Absorbed Charges, Boundary Conditions and 
Correlation Functions. 

A useful special case occurs when one has absorbed 
charges. These are charges that vanish into (or are in- 
jected from) absorbing nodes without altering the sys- 
tem dynamics. In mesoscopics for example, the absorb- 
ing nodes are metallic reservoirs. Formally, we divide 
the charges into those that are conserved and those that 
are absorbed: Q = {Q c , Q a }, where the subset of ab- 
sorbed charges Q a = {Q%} does not appear in H a p. We 
do the same for the corresponding counting variables: 
A = {A c ,A a }. Because H a p does not depend on Q a , 
these charges may be integrated out by integrating the 
action by parts, 

if dt' A a ■ Q Q = - i J dt' Q Q • A a 

+ *(A?.Q?-A?-Q?), (15) 

and then functionally integrating over Q" to obtain 
<5(A°), where S is a functional delta function. This im- 
mediately constrains the A a to be constants of motion 



» n— 1 „ 

E/(Q/,Qi,i)= dA J[ dQ k dA k exp 

J fe=i ^ 



n-l 

-i^k ■ (Qfc+i - Qfc) + AtH(Q k ,A k ) 

Lfe=o 

I 



6 



so the functional integration over A a becomes a normal 
integration, T>A a — ► dA a . The absorbed kinetic terms in 
the action may then be integrated to obtain 



U(Q f ,Q h t)= J dA a y"pQ c PA c exp{5(Q,A)} . 



(16a) 

S(Q,A)= / dt'[-iA c ■ Q c + (1/2) y^H a p{K-X fj )] 

J ° a/3 

-zA a • (Q a f - Qf). (16b) 

Often one is interested in the probability to transmit 
some amount of charge through each of the absorbing 
nodes. By applying a Fourier transform to Eq. I|16a|) 
with respect to Q a (t) — Q a (0) we remove the last term 
in Eq. (|16b|) and obtain the path integral representation 
for the characteristic function Z which generates current 
moments at every absorbing node 



Z(A a ) = ypQ c 2?A c exp{5(Q,A)}, 



(17a) 



5(Q,A) 



dt'[-iA c -Q c 



-(l/2)J2H aP (X a -Xp)].(17b) 



a/3 



Note that the counting variables A° enter the action 
(|17b|l only as a set of constant parameters. The initial 
condition in the path integral (|17|) is given by the initial 
charge states Q c (0). There is a choice of the final condi- 
tion: by fixing the final Q c (t) one obtains the distribution 
of the conserved charge subject to this constraint, while 
by fixing A c (t) the corresponding characteristic function 
is obtained. The choice of A c (t) = in Eq. (|17fl gives the 
characteristic function of the absorbed charge under the 
condition that the conserved charge is not being moni- 
tored, i.e. the final charge state is integrated over. There- 
fore \nZ becomes the generator of the FCS, defining the 
charge cumulants at the absorbing node, 

fin 1„ 7 

«[W)-Q*(o)] n »= ' - 



(18) 



A»=0 



d(lX%) 

In the long time limit, this quantity is proportional to 
time, independent of the details of the boundary condi- 
tions. 

Alternatively, in the short time limit one may calculate 
irreducible correlation functions of absorbed and con- 
served current fluctuations, I = Q. These correlation 
functions can be obtained by extending the time integral 
in (|3b|) to infinity, introducing sources^ in the action, 
S — > S + i J dt\(t) - , and applying functional deriva- 
tives with respect to \. Repeating the steps leading to 
Eqs. (|17fl . we find that variables A Q in the Hamiltonian 
in Eq. Ijl7b|l have to be shifted X a — » A Q -I- %a ■ Then, the 
irreducible current correlation function is given by 

S n lnZ\x] 



((i ai (h) •••/«„(*«)» = 



(19) 



With these correlation functions, one may calcu- 
late for example the frequency dependence of current 
cumulants. 34 



C. The Saddle Point Approximation. 

If the Hamiltonian has some dimcnsionless large pref- 
actor, then the path integral Q may be evaluated using 
the saddle point approximation, which is justified below. 
At the saddle point, (where the first variation of the ac- 
tion vanishes), we can write equations of motion analo- 
gous to the Hamiltonian equations of classical mechanics: 



tQ c = —H(Q c ,A), ik c 



dQ 



: ff(Q c ,A), (20) 



where H{Q C ,A) = (l/2)J2 a/3 H a [ 3 (Q c ;X a -Xp). There 
may be many saddle point solutions in general, and one 
has to sum over all of them. Eqs. (|2U|) are solved sub- 
ject to the temporal boundary conditions and generally 
describe the relaxation of the conserved charges from the 
initial state to a stationary state {Q c , A c } on a time scale 
given by tc, the dynamical time scale of the nodes. These 
stationary coordinates are functions of any external pa- 
rameters as well as the (constant) absorbed counting vari- 
ables A a . In the saddle point approximation, the action 
takes the form S = S sp + Sfi uc £ The term S sp is the 
contribution to the action from the solution of the equa- 
tions 120(1 . which describes the evolution of the system 
from the initial to the final state. The term S 



flue 



de- 



scribes fluctuations around the saddle point and is sup- 
pressed compared to the saddle-point contribution, if the 
Hamiltonian has a large prefactor (in analogy to the h- 
expansion of quantum mechanics). Physically, the valid- 
ity condition for the saddle point approximation is that 
there should be many (transporting) charge carriers in 
the nodes. For times longer than the charge relaxation 
time of the node, the dominant contribution is from the 
stationary state only, where the saddle-point part of the 
action is simply linear in time: 



S sp (Q,A)=tH(Q,A), i»r c 



(21) 



The linear time dependence of Eq. (|21|l indicates that the 
dynamics are Markovian on a long time scale. It is the 
fact that the contribution S sp emerges in a dominant way 
which makes the approach given here a powerful tool to 
analyze the counting statistics of transmitted charge. 

We now discuss the large parameter that justifies the 
saddle point approximation. The boundary conditions 
on the charge in the absorbing nodes fix a (dimension- 
less) charge scale of the system, 7. All charges in the 
network are scaled accordingly, Q — > 7Q. We make the 
assumption that there is a one parameter scaling of the 
Hamiltonian, H — > •yH. The time is also scaled by tc, 
the time scale of charge relaxation in the nodes. The di- 
mensionless action is now S — j dt'(—iQX + tcH). 
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The saddle point action is proportional to jt/rc, while 
the fluctuation contribution will be of order t/r c . We 
note that the parameter 7 is related to (though not nec- 
essarily the same as) the separation of time scales, tc/tq, 
needed to derive the path integral. For the mesoscopic 
conductors considered in the example section IV Bl of this 
paper, the charge scale is set by the maximum number of 
semiclassical states on the cavity involved in transport, 
7 = A^iNp 3> 1, the bias times the density of states 
at the Fermi level. On the other hand, for the chaotic 
cavity, t c /t = j/(G L + G R ), where G L ^ R > 1 arc the 
dimensionlcss conductances of the left and right point 
contact. 



D. Relation to the master Equation and Doi-Peliti 
Technique. 

The evolution operator J7(Q, Q', t) may be interpreted 
as a Green function of a differential equation which de- 
termines the propagation in time of an initial probability 
distribution T(Q). In the theory of stochastic processes, 
such a differential equation is called a master equation. A 
natural question that arises is the relationship of the for- 
malism presented here to other approaches to stochastic 
problems. 

The most general type of Markovian master equation 
for discrete states and discrete time is of the form 



r n (ife+l) — / Pnm(tk+l>tk)T m (tk) 



(22) 



where T m (tk) is the probability to be in state m at time 
tk and P nm is the transition probability from state m 
to state 71. The state is described by a vector n — 
(ni,...,njv) whose components are the charges n a of 
each node a. The Markovian assumption implies that 
tfe+i — tk — At is greater than the correlation time, tq. If 
we further assume that the probability to make a tran- 
sition to another state is small, P nrn <C 1 for n 5^ m, 
so that the transition probability is only linear in At, a 
transition rate W nm — P nm / At may be defined. It then 
follows that we may write a differential master equation, 



r«(t) = Y\w nm v m (t) - w m „r„(t)] 



(23) 



Eq. (|23|l is the starting point for the Doi-Peliti 
technique )2» where one formally maps the space of 
physical states to the Fock space of states \n) — 
{a\) ni . . . (<J N ) n *'\Q), where n is the number of charges. 
The entire state of the system is expressed by a vector 
1^) = Tln^ n \ n ) w hich weights the states \n) with their 
probabilities T n . Thus, the master equation may be 
interpreted as a many-body Schrodinger equation where 
the rates W mn are incorporated into a Hamiltonian in a 
second-quantized form. One may then write a coherent- 
state path integral over the variables a, and for this 



many-body quantum system and perform perturbation 
expansions along with the renormalization group^ This 
procedure eventually involves taking the continuum limit 
so the discrete charge states become continuous. 

Let us now consider how our formalism is related to 
the master equation or the Doi-Peliti technique. Accord- 
ing to the results of Sec. |HJ our stochastic path integral, 
Eq. © solves the continuum variable version of Eq. 11' I'D 
with the transition probabilities given by the one step 
propagator ?7(Q, Q', At). In general, the transition prob- 
abilities are neither small nor linear in time for At > tq. 
It is instructive nevertheless to consider the special case 
of processes where Htq <C 1, when we can expand the 
one-step propagator l|13|) to first order in At, 

U(Q, Q', At) pa <5(Q-Q') 

+ At J dAe~ lA -( Q - Q '>H(Q',A). (24) 

Defining the Fourier transform of the generating func- 
tion as H(Q, Q'), the differential equation governing the 
evolution of a probability distribution of charges T(Q) is 
then 



f(Q,t)= / dQ'ff(Q,Q')r(Q',t). 



(25) 



Comparison with the continuous version of the master 
equation l|23|) . 

f(Q,t) = J dQ' [W(Q,Q')nQ',t)-W(Q',Q)T(Q,t)}, 

(26) 

indicates that H is related to W. The Hamiltonian may 
be expressed in terms of the transition kernel^ as, 



i?(Q',A) = / dQ 



i(Q-Q') A _ 1 



W(Q,Q'), (27) 



where the normalization of probability is expressed by 
H(Q',0) = 0. Eq. I)27|) is an important result, because 
it allows the conversion of the master equation (|26|) into 
the stochastic path integral 10). 

We would like to stress that our formalism is not simply 
equivalent to the differential master equation (|26|l (and 
therefore the Doi-Peliti technique), but that it allows the 
treatment of a complementary class of problems. Our for- 
malism assumes effectively continuous charge, and thus 
cannot resolve effects due to the discreteness of charge on 
the nodes. Such effects are present in the master equa- 
tion (|23|) . In contrast, the differential master equation 
assumption, Htq <C 1 (which simply states that transi- 
tion probabilities are small in the time interval To) is not 
required. Our formalism is especially important when 
this is not the case, i.e. Htq ~ 1. 

This is illustrated by the simple example from meso- 
scopics of two metallic reservoirs connected by a sin- 
gle electron barrier with hopping probability p and bias 
A/x at zero temperature. For a time interval At larger 
than the correlation time tq = ft/A/j (the time scale for 
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for this arrangement is 



R 



FIG. 2: A one dimensional lattice of nodes connected on both 
ends to absorbing reservoirs. This situation could represent 
a series of mesoscopic chaotic cavities connected by quantum 
point contacts. 



an electron wavepacket to transverse the barrier), At/ tq 
electrons approach the barrier and either are transmitted 
or reflected. Mathematically, this is a classical binomial 
process with the generator 



5=(Ai/ro)ln[l+p(e leA -I)]. 



(28) 



As this action is the starting point of many mesoscopic 
implementations of the formalism, it is an important ex- 
ample. Since the action is proportional to the large pa- 
rameter At/To > 1, for p ~ 1 the expansion of exp(5) 
to first order in At is strictly forbidden, effectively not 
allowing a first order differential master equation. Only 
in the limit p <C 1, (i.e. when Eq. (|28|l describes a Pois- 
sonian process) may the logarithm be expanded to first 
order. This suggests that Eq. I)26|l describes the slow 
dynamics of systems whose fast transitions are Poisso- 
nian in nature. A more general type of dynamics such 
as the binomial distribution may only be found using the 
continuous charge state master equation in discrete time 



III. THE FIELD THEORY 

From the stochastic network, Fig. ^ it is straightfor- 
ward to go to spatially continuous systems as the spacing 
between the nodes is taken to zero. The goal is to intro- 
duce a Hamiltonian functional h(p, A) whose arguments 
are the charge density p and the counting field functions 
A, that are themselves functions of space and time. We 
may then replace (1/2) Y^, a g — > / dzh(p.X). Our 
description is local, so in the model each node is only 
connected to its nearest neighbors. We first derive the 
one dimensional field theory with one charge species in 
detail, and then generalize to multiple dimensions and 
charge species. 

Consider a series of identical, equidistant nodes sepa- 
rated by a distance Az. This nodal chain could repre- 
sent a chain of chaotic cavities, Fig. |3 in a mesoscopic 
context^a The sum over a and (3 becomes a sum over 
each node in space connected to its neighbors. The action 



Jo 



A Q Q Q + H(Q a , Q„_i; A Q — A Q _i)} , 

(29) 

where for simplicity we have chosen real counting vari- 
ables, iX a — > A Q . The imaginary counting variables 
will be restored at the end of the section. The only 
constraint made on H is that probability is conserved, 
H(X a — A„_i) = for A Q = A a _i. We now derive a lat- 
tice field theory by formally expanding H in X a — A Q _i 
and Q a — Q a ~i- Only differences of the counting vari- 
ables will appear in the series expansion, while we must 
keep the full Q dependence of the Hamiltonian. If there 
are N ^> 1 nodes in the lattice, for fixed boundary condi- 
tions the difference between adjacent variables, A Q — A Q _i 
and Q a — Q a -i will be of order 1/N, and therefore pro- 
vides a good expansion parameter. The expansion of the 
Hamiltonian l|29|) to second order in the difference vari- 
ables gives 



H = 



OH 1 d 2 H 2 

^-(A Q -A Q _ 1 ) + -^— (A Q -A Q _x) 



d 2 H 

dQ a dX a 



(Qa 



la-l 



)(A Q — A Q _i) . 



(30) 



where the expansion coefficients are evaluated at X a — 
A Q _i and Q a — Q a -i and are functions of Q a -i- Terms 
involving only differences of Q a — Q a -i are zero be- 
cause H(X a — A Q _i) = for A a = A Q _i. All terms in 
Eq. I|3UI) need explanation. First, the expression dH/dX a 
is the local current at zero bias (because the charges 
in adjacent nodes are equal) which will usually be zero. 
There may be circumstances where this term should be 
keptj2£ but we do not consider them here. The term 
d 2 H/dQ a dX a = —G(Q a -i) is the linear response of the 
current to a charge difference. Hence, G is the gener- 
alized conductance^ of the connector between nodes a 
and a — 1. d 2 H/dX a — C(Q a ~\) is the current noise 
through the same connector because H is the generator 
of current cumulants. 

We are now in a position to take the continuum limit 
by replacing the node index a with a coordinate z, intro- 
ducing the fields Q(z), X(z), and making the expansions 

A Q - A Q _i -> X'Az + (l/2)A"(Az) 2 + 0(Azf, (31a) 
Q a - Q a -i -^Q'Az + (l/2)Q"(Az) 2 + 0{Azf. (31b) 

The action may now be written in terms of intensive fields 
by scaling away Az, 



H — > h(p, A)Az, Q a — > p(z)Az, 
G Q (Az) 2 -> D(p), C a Az^F(p). 



(32) 



and taking the limit ^2 a H — > J dzh(p,X). One may 
check that expanding the Hamiltonian to higher than sec- 
ond order in Az will result in terms suppressed by powers 
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of Az/L and consequently vanish as Az — > 0. This scal- 
ing argument for the field theory is analogous to Van 
Kampen's size expansion^ Though the lattice spacing 
Az does not appear in the continuum limit, it provides 
a physical cut-off for any ultra-violet divergences that 
might appear in a loop expansion. 

These considerations leave the one dimensional action 

as 



time, thus obtaining 



S 



dt' / dz 

n JO 



\p + Dp>\'-~F(\'f 



(33) 



Here D is the local diffusion constant and F is the lo- 
cal noise density which are discussed in detail below, ft 
is very important that these two functionals D, F are 
all that is needed to calculate current statistics. Classi- 
cal field equations may be obtained by taking functional 
derivatives of the action with respect to the charge and 
counting fields: SS/Sp(z) — SS/SX(z) = to obtain the 
equations of motion, 



DA 



p = [-FX' + Dp'} 1 . (34) 



From the charge equation, one can see immediately that 
the term inside the derivative may be interpreted as a 
current density so that local charge conservation is guar- 
anteed. We have to solve these coupled differential equa- 
tions subject to the boundary conditions 



P (t,0) =p L (t), p(t,L)=p R (t), 
A(t,0) = X L (t), X(t,L) = X R (t), 



(35) 



where PLif), Pn(t), Ai(i), and X R (t) are arbitrary time 
dependent functions. Functions px,(i) and p R (t) are the 
charge densities at the far left and right end of the system 
which may be externally controlled. Functions Al(£) and 
X R (t) are the counting variables of the absorbed charges 
at the far left and right end which count the current that 
passes them. 

Once Eqs. (|34|l are solved subject to the bound- 
ary conditions l|35|l . the solutions p(z,t) and X(z,t) 
should be substituted back into the action (|33|) and 
integrated over time and space. The resulting func- 
tion, S a p[pi(t),p R (t),\ L (t),\ R (t),t,L] is the generating 
function for time-dependent cumulants of the current 
distribution. Often, the relevant experimental quan- 
tities are the stationary cumulants. These are given 
by neglecting the time dependence, finding static solu- 
tions, p = A = Q, and imposing static boundary condi- 
tions. Similarly to section III Dl we can also introduce 
sources J dtdz x(z,t)p(z,t) and calculate density correla- 
tion functions. 

To estimate the contribution of the fluctuations to the 
action, it is useful to define dimcnsionlcss variables. The 
boundary conditions pl, and p R provide the charge den- 
sity scale po in the problem, so we define p(z) = pof(z), 
where / ~ f is an occupation. We furthermore rescale 
2 — > Lz, and t — * Tot, where td = L 2 /D is the diffusion 



S 



-Lp 



dt' 

^0 



dz' 



Xf + f'X' 



F 



2Dp Q 



(X'f 



(36) 

We assume that the combination F/Dpo is of order 1. 
From Eq. 1|36[) . the dimensionless large parameter is 7 = 
poL 3> f , i.e. the number of transporting charge carriers. 
As in Sec. Ill CI the saddle point contribution is of order 
jt/TD, while the fluctuation contribution is of order t/rjj. 

Repeating this derivation in multiple dimensions with 
iV" charge species p = {pi(r)} and counting fields A = 
{Aj(r)}, i — 1, . . . , N yields the action 

S = - [ dt' [ dr [Ap + VADVp- (1/2)VAFVA], 
Jo Jn 

(37) 

where tensor notation is used and we have introduced 
Fij = d\idxjh and Dij = —dpfix^h as general matrix 
functionals of the field vector p and coordinate r which 
should be interpreted as noise and diffusion matrices. If 
the medium is isotropic, then the vector gradients simply 
form a dot product. It should be emphasized that the 
vectors appearing are vectors of different species of charge 
fields, as all node delimitation has been accounted for 
in the spatial integration. The functional integral now 
runs over all field configurations that obey the imposed 
boundary conditions at the surface dVL. Classical field 
equations may be formally obtained by taking functional 
derivatives of the action with respect to the charge and 
counting fields as in the ID case. 

As in any field theory, symmetries of the action play 
an important role because they lead to conserved quan- 
tities. We first note that the Hamiltonian h(p, \7p, VA) 
is a functional of VA alone with no A dependence. This 
symmetry is analogous to gauge invariance, and leads to 
the equation of motion 



P + V-j = 0, j = -D\7p + FVA, 



(38) 



which can be interpreted as conservation of the condi- 
tional current j. The next symmetry is related to the 
invariance under a shift in the space and time coordi- 
nates {Sr 7 St}. This symmetry leads to equations analo- 
gous to the conservation of the local energy/momentum 
tensor. 4 S We do not explicitly give this quantity because 
it is rather cumbersome in the general case. However, for 
the stationary limit (where p and A vanish) and for sym- 
metric diffusion and noise tensors, the one charge species 
conservation law is relatively simple and is given by 

VmTmn = , (39a) 

i n 

T mn = im(V„A) - (V n p) (DVX) m ~ h S mn . (39b) 

For the special case of a one dimensional geometry, 
the Hamiltonian itself is the conserved quantity (see 
Sec.ES). 
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In the continuum limit, all terms of higher order in A 
are suppressed so that the action is quadratic in the A 
variables. This fact may be viewed as a consequence of 
the central limit theorem and confirms the observation 
made by Nagaev that local noise in the mesoscopic dif- 
fusive wire (see Sec. IV A} is Gaussian; 1 - To further clar- 
ify the physical meaning of D and F, and also to make 
connection with previous workjS we restore the complex 
variables, A — > iA, and make a Hubbard-Stratronovich 
transformation by introducing an auxiliary vector field 

exp{-(l/2)VAFVA} 

= (dctPyi J Vv exp{-{l/2)vF~ 1 v + ivVA}. (40) 

We may then integrate out the A variables, taking ac- 
count of the boundary terms to obtain, 

U = cx p{^ dt'J ds-{iA a 3)\Jv P VvS{p + V -J) 

x (detF)-^expj-i^ dt'J^dr'vF' 1 ^ , (41) 

where the 6 above is a functional delta function, imposing 
the Langevin equation 

p + V • J = , (42a) 
J = -DVp + v, (42b) 

with a current noise source v, whose correlator— is given 
by 

(i/(r, t)v(r', t')) = S(t - t')8(r - r')F(p). (42c) 

J may be interpreted as the physical current density [not 
to be confused with the conditional current density (j2BJ)] 
so that local current conservation is guaranteed, and the 
(detF) -1 / 2 serves to normalize the v probability distri- 
bution. The role of the boundary term is to count the 
current J flowing out of the boundary with the count- 
ing variable A a , which serves as a Lagrange multiplier. 
This formula gives an immediate translation between the 
Langevin approach and full counting statistics, a connec- 
tion not previously known. The algorithm is as follows: 

1. Given a Langevin equation of the form (|42|l . write 
the average of the boundary term with source A° as 
a path integral (|41|l over noise and density fields— 

2. Introduce an auxiliary field A that takes on the 
value A a at the boundaries and represents the delta 
function in Eq. (|41|l imposing current conservation 
(|42a|) in Fourier form. 

3. Integrate out the Gaussian noise to obtain an action 
of the form of Eq. (|37|) . 

4. Find where the first variation of the action is zero 
and solve the equations of motion subject to the 
boundary conditions. 



5. Insert the solutions back into the action, and do the 
space and time integrals. The answer is the current 
cumulant generating function. 



IV. PERTURBATION THEORY 

We have shown in Sec. Ill Gl that a large number of par- 
ticipating elementary charges justifies the saddle point 
approximation for the generator of counting statistics. 
While the generator may sometimes be found in closed 
formji in general, it has no compact expression and the 
cumulants should be found separately at every order. 
This may be done by expanding S sp (Q 1 A,x) as a series 
in x and solving the saddle point equations to a given or- 
der in x directly. However, there is another approach for 
evaluating the higher cumulants, the cascade diagram- 
matics representing higher-order cumulants in terms of 
the lower ones. It has been introduced by Nagaev in 
the context of mesoscopic charge statistics in the diffu- 
sive wirei^ and later extended to the chaotic cavity, 20 
but without proof. The basic idea is that lower order 
cumulants mix in to yield corrections to the bare fluctu- 
ations of higher order cumulants. This method was used 
successfully in Ref. |43j to explain the recent experiment 
of Ref. 44. In this section, we demonstrate that these 
rules follow naturally from the stochastic path integral 
in the same way as Feynman diagrams follow from the 
quantum mechanical functional integral. In Sec. II V CI we 
present another (simpler) method for computing cumu- 
lants based completely on differential operators obtained 
from the Hamiltonian equations of motion. In Sec. IIV Dl 
we generalize the cascade diagrammatics to an arbitrary 
network, and to the case of time-dependent correlators. 



A. The Principle of Minimal Correlations. 

To motivate the cascade diagrammatics, we refer to a 
specific physical system (see the inset of Fig.[SJ), the meso- 
scopic chaotic cavityi For the purposes of this section, 
the cavity is a conserving node carrying charge Q, the 
electronic reservoirs correspond to the left and right are 
absorbing nodes, and the two point contacts are the con- 
nectors described by Hamiltonians Hl,Hr (see Fig. [3J. 
Although a detailed description of this system is given in 
Sec. IV Bl we would like to mention that the mesoscopic 
cavity is described by an electron distribution function 
/, which is fluctuating around its mean value, /o- The 
actual electrical charge in the cavity Q and the occu- 
pation / are related via the large parameter 7 through 
Q = 7(/ — Jo), where 7 = ApNp 3> 1 (the density 
of states at the Fermi energy Np times the bias A/i) is 
the maximum possible number of electrons on the cavity 
which contribute to the transport (see Sec. Ill C|) . 

The cascade approach builds on the principle of mini- 
mal correlations developed in Ref.^3 The point contacts 
create bare noise ((/£)) = d 2 H L / '(OiXl) 2 , and {{Ir}} = 
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contains a contribution from fluctuations of the charge 
in the cavity that couples back into the current fluctua- 
tions. The factor of 3 comes from the fact that there are 
3 independent currents that the charge fluctuation may 
be correlated with. For higher cumulants, there will be 
more cascade corrections that may be represented in a 
diagrammatic form^SS 



FIG. 3: Network representing a chaotic cavity. The state of 
the internal node is described by the variable Q, the charge on 
the cavity. The statistics of the connectors are characterized 
by the two generating functions Hl,r- 



d 2 H R / (di\n) 2 with no correlation, ((IlIr)) = [see Eq. 
|JSJ]. However, for times longer than the average dwell 
time of electrons in the cavity, the current conservation 
requirement imposes "minimal correlations" on the fluc- 
tuations of the physical currents II and Ir, which can 
be expressed in the form of the Langevin equations, 



Il = Il- G l Q, 



Ir = h. 



GrQ, 



(43) 



where Il,r are now the sources of bare noise, Gl,r are 
the generalized conductances of the left and right point 
contact, and Q is the fluctuating charge in the cavity. 
Current conservation of the physical currents, II = Ir. = 
/, can now be used to obtain 



I = 



GrJl + GlIr 
Gl + Gr 



Q = 



Il-Ir 
Gl + Gr 



(44) 



Combining powers of / and Q and averaging over the 
bare noise, we obtain the minimal correlation result 
for arbitrary cumulants ((Q k I l )) m - In particular, using 
{(IlIr)) = 0, we find the second cumulant of current 

:„17.18 



<(^}} = ((I 2 ))r 



Gr((II))+GI((I r )} 

(G l + Gr) 2 



(45) 



where the subscript m denotes the minimal correlation 
result. We stress that the bare correlators ((/£ R )) are 
fully determined by the average occupation function / 
of the cavity. 

This example demonstrates that a simple redefinition 
of the current fluctuations makes it straightforward to 
find the noise. Therefore, it came as a surprise^ that the 
minimal correlation approach is not sufficient to correctly 
obtain higher-order cumulants of current. The reason 
for the failure of the minimal correlation approach has 
been found recently by Nagaev^ who showed that from 
the third order cumulant on, there are "cascade correc- 
tions" to the minimal correlation result, which may be 
interpreted as "noise of noise" . For example, the third 
cumulant of current through the mesoscopic cavity^ 



«J 3 » = «J 3 )) ro + 3«/Q)) r 



d_ 

dQ 



(46) 



B. Derivation of Diagrammatic Rules. 

We now present a derivation of these diagrammatic 
rules for a single node attached between two absorbing 
nodes. Generalizations to an arbitrary network will sub- 
sequently be given in Sec. IIVDI As we have shown in 
Sec. Ill CI the charge scale imposed by the boundary con- 
ditions, 7, gives a dimensionless large parameter which 
justifies the saddle point approximation of the path inte- 
gral, so that fluctuations around the saddle point are sup- 
pressed by I/7. In the diagrammatic language, we will 
show that loop diagrams are suppressed by the same fac- 
tor I/7. The diagrammatic approach given here is based 
on perturbation theory originally developed in quantum 
mechanics ~ 

Consider the path integral expression of the generating 
function for the charge absorbed in the left (L) and right 
(R) node: 

Z(xl,Xr) = I VQV^ exp( / dt'[-iQ\ 



+H(Q,X, X l,Xr)}}, (47) 



where H = H L (Q, X - xl) + H R (Q, X R ~ A). The per- 
turbation theory is formulated as follows. First, the ex- 
ternal counting variables are set to zero, xl = Xr = 0. 
The Hamiltonian H — > Hl(Q, X) + Hr(Q, — A) has a sta- 
tionary saddle point located at {Qo,\o} that we wish 
to define as the origin of coordinates. The probability 
distributions of transferred charge are normalized, so 



d%H L , R (Q,\)\ x = 



0, 



Vn. 



(48) 



In particular, d Q H L (X)\\=o = dQH R (X)\ x =o = 0, and 
therefore A = 0. Next, d lX H(X)\ x=0 = (I L (Q)) - 
(Ir(Q)) = 0, since Hl and Hr are the generators of 
the left and right current respectively. Therefore, Qo is 
fixed as the charge in the node such that left and right 
connector currents are equal on average. The stability of 
the saddle point is guaranteed by the fact that the bare 
noise correlators, ((Ilr)}, are positive. The derivatives 
di\dQH L = -G L , d l xd Q H R = -Gr define the gener- 
alized conductance of each connector, where the current 
flows from left to right in both connectors. 

The principle of minimal correlation plays an impor- 
tant role in the cascade diagrammatics. We will show 
that this principle is equivalent to exploiting certain free- 
doms in the path integral in order to postpone the cas- 
cade corrections to third and higher order cumulants. In 



12 



the long-time limit, t 3> tq (where 1/rc = Gl+Gr is the 
relaxation rate of the charge in the node), the absorbed 
current is conserved, Ir — II- Therefore, the current 
through the node can be defined as weighted average of 
the left and right connector currents I = (1 — v)Il + vIr, 
where v is arbitrary constant. The corresponding count- 
ing variable x is introduced by substituting xr — vx and 
Xl = (v — l)x- Consider now the second derivative 



d 2 H 



dixdQ 



(v-l)G L + vGj 



(49) 



We may set it to zero by fixing v — Gl/{Gl + Gr). 
This is equivalent to imposing conservation of current 
fluctuations as in Eq. (|44ll . If we consider further the 
derivative 



d 2 H 



di\dQ 



— —(Gl + Gr), 



(50) 



we have the freedom to scale A to make the right hand 
side of Eq. I|5U|I equal to —1 [this scaling only alters the 
X independent prefactor of Eq. I|47(l ]. The Hamiltonian 
takes the new form 



H = Hl Q, 



GrX + ^ 
Gl + Gr 



Hr[Q, 



Gl + Gr 



(51) 



We refer to these new variables as minimal correlation 
coordinates and will see that they simplify the diagram- 
matic expansion. 

Define 5Q(t) = Q(t) - Q and 5X(t) = X(t) - A . If 
we expand the Hamiltonian in a power series in x, 5Q, 
and <5A, the terms linear in SQ and <5A vanish at the sad- 
dle point, as well as the (5Q) 2 coefficient by Eq. H48|) 
with n — 2. As argued above, in the minimal correlation 
coordinates, di\dqH{Q 0l Ao) = — 1. With these transfor- 
mations, we may split the action S as 

S = S + f dt'V{t'), So = -i I dt' SX(t c SQ + SQ), 



(52) 

where V represents the rest of the H power series and will 
be treated perturbatively. It should be emphasized that 
V is a general nonlinear function of SX, so unlike most 
quantum examples, the full momentum dependence must 
be kept. 

In order to formulate the perturbation theory, we add 
two sources, J and K to the action, S — > S + J dt'[JSQ + 
iKSX], so that any average of a function of the variables 
SQ, SX may be evaluated by taking functional derivatives 
with respect to the sources J, and K, and then setting the 
sources to zero. In particular, for the generating function 
we can write 



Z(x) = JvQVXexp^ dt'V(SQ,5X,x)^exp^S + dt' [ J5Q + iKSX] j 
= exp|^ dt'V (^>^>x) jy^Q^Aexp js +^ dt'[J8Q + iK8X]^ 



J,K=0 



J,K=0 



(53) 



Using Sq from Eq. (|52(l we evaluate the integral over Q 
and A and obtain: 



where the functional W(J, K) is 



j,k=o 
(54) 



W(J, K) = cxp 



dt'dt"J{t')D{t',t")K{t") 



(55) 

The operator D = (rcdt + 1) 1 is the retarded propaga- 
tor, and may be found explicitly by inverting the kernel 
in frequency space, 



D(t,t') 



du: e^'-*') 
-oo 2tt -ircui + 1 

= Tc'e^-tOexpKi-O/rcL (56) 



It describes the relaxation of the charge Q(t) to the sta- 
tionary state Qq with the rate 1/rc = Gl + Gr. 



Expanding the exponential in Eq. I|54|l and taking the 
t 3> tq limit, we arrive at the following expression for 
the nth current cumulant 



<(/"}} = r 



S" 



oo 

y — 

S(iY) n ^— ' m! 



dt ' V {5J>5lK^ 



xW(J,K) 



X = .1=K = 

connected 



(57) 



According to the linked cluster expansion^ by consid- 
ering lnZ(x) rather than Z{x)i we have eliminated all 
disconnected terms. In order to compare with the results 
of Ref. 1201 we introduce a new notation by defining 



dU(Q k I l )) m ee 9^ a Cf(Qo,A ,x = 0). 



(58) 



Here ((Q k I l )) m is the irreducible correlator expressed in 
terms of the noise sources, i.e. the minimal correlation 
cumulant. In this notation, the expansion of V in a Tay- 
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■--> = 1, • > =0. 

(c) (d) 

FIG. 4: (a) An n-point current cumulant. (b) The vertex 
connecting I external lines with j internal Q lines and k in- 
ternal A lines, (c) The propagator connecting A to Q, equal 
to 1 in the stationary limit, (d) The vanishing vertex 8q{I) 
in minimal correlation coordinates. 

lor series of all variables takes the form: 
V(SQ,SX, X ) = ^-L-<P Q ((Q k I l )) m 

x[5Q(t)}1[i5\(t)} k [i X } 1 . (59) 

Inserting the expansion Eq. (|59|l into the formula for the 
current cumulants Eq. (|57f) gives the formal solution to 
the problem. From the form of W(J, K) and V, we can 
immediately read off the diagrammatic rules with the in- 
ternal lines given by the propagators (|56[) . and the expan- 
sion coefficients d J Q((Q k I l )) m playing the role of vertices. 

The following simplifications can be done before the 
rules are finally formulated. First, it is straightforward 
to see that loop diagrams are suppressed by powers of 
7 . Indeed, according to our single-parameter scaling 
assumption, the action (|52|l has a large prefactor 7, which 
can be explicitly displayed, S — > 7 S, by rescaling the 
charge, Q — » 7 Q. Then it becomes clear that each prop- 
agator D, represented by an internal line, comes with a 
factor of 7 _1 . Each vertex comes from V and therefore 
has a factor of 7. If a diagram has I internal lines, E 
external legs, V vertices and L loops, it will come with 
a total 7 power of V — I. Furthermore, Euler's formula 
tells us that V + L — I = 1. Therefore, diagrams with no 
loops ("tree" diagrams) come with a power of 7, while 
loop diagrams are suppressed by the number of loops, 
7 1_L . From now on we will concentrate on tree-level di- 
agrams, since they represent current cumulants at the 
level of the saddle-point approximation. 

Second, in the long time limit, t 3> tc, each propagator 
(|56|l integrated over time gives 1. As a result, since every 
vertex is connected to at least one other vertex, all the 
time integrals together simply give a factor of t, and the 
time dependence cancels on the right hand side of the 
Eq. (|57|l . There are no time integrals in the vertices and 
the propagators just give a factor of 1 as in Ref. 12(3 . 
We are now able to formulate the diagrammatic rules for 
high-order current cumulants: 



1. The nth order cumulant ((/")) is a connected n- 
point function of n external legs / represented by 
solid arrows (see Fig. 0^). 

2. The external legs must be connected by using ver- 
tices (see Pig. St>) and linking internal dashed lines 
to internal dashed arrows. 

3. The vertices d-Q{{I l Q k )) m are represented by a cir- 
cle with I external legs, k internal outgoing dashed 
lines, and j internal incoming dashed arrows (see 
Fig.B>). 

4. Multiply each diagram by the number of inequiva- 
lent permutations (NIP). 

Formally, the vertices d 3 Q((I l Q k )) m are the expansion co- 
efficients in (|59p. However, it is important to note that 
they can also be easily evaluated by solving the Langevin 
equations l|43ll and expressing the minimal correlation 
cumulants ((I l Q k )) m in terms of cumulants of the noise 
sources, ((1^)) and ((Ir ))• Some vertices are zero, 
d p H/dQ p (Qo, Ao)| x =o = because of probability conser- 
vation, but other may or may not be zero depending on 
the physical system. Here, the advantage of the min- 
imal correlation coordinates is made clear: the vertex 
9q((/)) ot = 0, and therefore any diagram that contains 
this vertex is zero (see Fig. QJi). 

To obtain the overall prefactor of a diagram, one can 
write out all the numerical constants and count the num- 
ber of different ways of producing the same diagram, ? 2 
For example, there is the n\ from the \ derivatives, the 
1/m! from the Taylor series of e , a binomial coefficient 
from expanding V m , and the 1 / ( j ! A; ! Z ! ) from every vertex 
with j + k + l attachments for the different lines. To com- 
pensate these factors, we have to do the combinatorics of 
the number of equivalent terms: interchange the vertices, 
find the number of different placements of lines on a ver- 
tex, etc. Often, the number of permutations of the n 
external legs will cancel the m!, and the j\k\l\ number of 
permutations of the internal legs attaching to the vertex 
will cancel that factor arising from the Taylor expansion. 

Rather than making this expansion, there is a sim- 
pler method which exploits these cancellations given by 
counting the number of inequivalent permutations of the 
diagram (NIP). The NIP of the diagram is defined by how 
many ways the external legs of the diagram may be rela- 
beled, such that the diagram is not topologically equiv- 
alent under deformation of the external legs. In other 
words, a diagram with n external legs has n\ ways of la- 
beling them. If this diagram with a given labeling of the 
legs may be topologically deformed to give the diagram 
back with a different labeling, these two sets of labclings 
are equivalent permutations. If we write out all the dif- 
ferent labelings the external legs can have, and cross out 
every labeling that is an equivalent permutation of an- 
other, then the number of labelings that remain is the 
NIP. This number is most easily found by dividing n! by 
the number of equivalent permutations of the diagram. 
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(a) (b) (c) 

FIG. 5: Tree level contributions to the third cumulant of 
transmitted current. 



The number of equivalent permutations of the diagram 
is also called the symmetry factor of the diagram. 

We illustrate these two approaches with the third cu- 
mulant, see Fig. |SJ With the simplifications discussed 
above, these diagrams may be written as 

((I 3 )) = ((I 3 )) m + 3«JQ)) m ^((J 2 )) m 

+ 3«/Q»^^«/» m . (60) 

Note that diagram (c) does not appear in Ref. 113, be- 
cause it happens to vanish for the chaotic cavity [see also 
Eq. H46|) ]. Referring to the formula (|57|) . the contribu- 
tions in Eq. I|60() are from m = 1,2,3 respectively. Each 
diagram must have a x 3 term in the expansion. We first 
show the combinatorial method to obtain the prefactor: 
Diagram (a) has a factor of 1/3! from the number of 
permutations of the x variables, canceling the 3! from 
the x derivatives. Diagram (b) has a factor of 1/2! from 
the number of permutations of the \ variables, a factor 
of 1/2! from the Taylor series of the exponential, a factor 
of 2 from the binomial expansion of V 2 , and the 3! from 
the x derivatives, leaving a factor of 3. Diagram (c) has 
a factor of 1 /3! from the Taylor series of the exponential, 
a factor of 3 from the binomial expansion of V 3 , a factor 
of 1/2! from the number of permutations of the 5Q vari- 
ables, a factor of 2 from the functional derivatives acting 
on W, and the 3! from the x derivatives, leaving a factor 
of 3. The NIP is simpler to derive: We divide the number 
of permutations of the external legs, m!, by the number 
of equivalent permutation of the elements of the diagram 
that leave it unchanged. The number of equivalent per- 
mutations of diagrams (a,b,c) are 3!, 2!, 2!, leaving the 
overall factors 1, 3, 3. 




FIG. 6: Three examples of diagrams contributing to the 
fourth cumulant. 



The computation of these diagrammatic contributions 
is best understood by a little practice on some examples. 
Consider three of the diagrams that contribute to the 
fourth cumulant drawn in Fig. [f)| The diagrams symbol- 
ically represents the combinations: 

(«) = ^((I 2 ))m^((Q)) m ((IQ)) 2 nl (61a) 
8 3 

(b) = QQ§((T)U((IQ)) 3 m , (61b) 

(C) = <(Q 2 » m (^2«'»m) {{IQ)?m- (61C) 

To figure out the numerical prefactors, we divide 4! (4 is 
the number of external legs) by the symmetry factor of 
the diagram. We first consider the symmetry factor of 
(a): The upper two legs may be flipped, and the lower 
two legs may be independently flipped where the dotted 
arrows join without altering the topology of the diagram. 
Therefore, the symmetry factor is 2 x 2 = 4, and the 
NIP is 4!/4 = 6. Moving on to diagram (b), the three 
lower legs may be permuted amongst themselves to give 
a symmetry factor 3!, and therefore the NIP is 4!/3! = 4. 
Finally, diagram (c) may be flipped about its center for 
a symmetry factor of 2, giving a NIP of 4!/2 = 12. 

C. Operator approach. 

In the stationary limit, t 3> Tc, the action takes the 
form S = tH(Q, A, x) so that the evaluation of the cumu- 
lant generating function reduces to finding the station- 
ary point of the Hamiltonian H as a function of vari- 
ables A and Q. This can be done by solving equations 
8qH = and d\H = 0. The generating function is then 
obtained by substituting the solutions {Q, A} into the 
Hamiltonian. In the previous section we have shown that 
this problem can be solved using path integral methods, 
and the solution can be represented diagrammatically. 
In the next section we will exploit the full strength of 
the path integral formalism in order to generalize the di- 
agrammatics to an arbitrary network, and for the case 
of time-dependent charges. However, in the stationary 
limit, the conceptual simplicity of the problem of find- 
ing the stationary point of the function H indicates that 
there should exist a simple iterative procedure for evalu- 
ating the cumulants up to a given order. In this section 
we use classical mechanics methods to prove that this is 
indeed the case. 

We first make the variable transformation iX — > A, and 
i\ —> Xi so that the Hamiltonian becomes a real func- 
tion. For x = the saddle point is located at {Qo,Xo}- 
For non-zero x the saddle point moves to a new posi- 
tion {Q,A}, which depends on \, and the Hamiltonian 
H(Q,X,x) becomes the generator of cumulants of the 
current, 

((/")) = d"ff(Q, A, X )/dx n \x=o- (62) 
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By expressing the total x derivative in terms of partial 
derivatives, the average current can be written as 

(I) = (d x + Q'dq + \'d x )H{Q, A, x)\ {x =o,Q M}, (63) 

where Q' = dQ/dx,X' = dX/dx are x dependent. We 
wish to eliminate the functions Q' and A' and to express 
the cumulant in terms of the partial derivatives of H. 
This is done by applying a total derivative to the equa- 
tions of motion: [dgH]' = [d\H]' = and leads to two 
equations for Q' and A' which may be solved, 



Q' = 



{d x H,d x H} 

{d Q H,d x Hy 



X' = - 



{d Q H,o x H} 

{d Q H,d x HY 



(64) 



where {A, B} is the Poisson bracket, defined as {A, B} = 
OxAOqB — OqAOxB. The solutions have to be inserted 
into the Eq. (|BH) . 

The advantage of this representation is clear: Now the 
right hand side of the Eq. (|63fl (before taking the x = 
saddle point) depends only on variables A, Q, and X- 
Therefore, we can apply the procedure again in order 
to express the high-order cumulant in terms of partial 
derivatives. This procedure solves the problem by giving 
a single operator, 



D = <9 V 



{d x H,d x H}d Q -{d Q H,d x H}d, 
{d Q H,d x H} 



(65) 



which, being applied n times to a given Hamiltonian H 
and evaluating the resulting expression at the x = sad- 
dle point, gives cumulants of current: 



((/")) = D n H(Q, X, x)| {x =o,Q ,Ao}- 



(66) 



This approach is obviously more simple compared to the 
diagrammatic method, since in the diagrammatics, af- 
ter drawing all of the diagrams, they have to be eval- 
uated individually by taking many partial derivatives of 
the Hamiltonian and evaluating them at the x — saddle 
point. With this new approach, given the Hamiltonian 
H, the operator D may be constructed (|65|l and with a 
mathematical program, an arbitrary cumulant may be 
easily computed 

It is easy to see the importance of the minimal correla- 
tion coordinates in this solution. After applying D several 
times, the derivative quotient rule generates a large num- 
ber of denominators, {dQH,d\H} = (dQd\H)(d\dQH) — 
(dQdQH)(d\d\H). At x = 0, as we argued previously, 
dqdqH = 0, and it is possible to change coordinates so 
that dqd\H = — 1. As a result, the denominator in i|6T)|) 
is equal to 1, which greatly simplifies the expansion. Fi- 
nally, we would like to stress that the operator approach, 
introduced in this section for the one node case, can be 
easily generalized to a network. 



D. 



Network Cascade Diagrammatics: Correlation 
Functions. 



Consider now a general network. In the Sec. lIVBI we 
saw that the dominant contribution to Eq. I|47|l arises 



from tree-level diagrams. On time scales t > tc, the 
time dependence drops out, and the current cumulants 
are static. We now generalize the diagrammatic rules 
presented in the section IIVBI to investigate time- and 
node-dependent correlation functions of conserved and 
absorbed charges, Eq. Q19[l. To define the network, we 
must arbitrarily label the current flow, yielding a directed 
network. By doing so we fix the signs of the elements 
H a f3 = —Hp a of the Hamiltonian. In particular, the 
elements of the generalized conductance matrix G, 



G 



d 2 H 



af3 



d{iX a )dQ{j 



(67) 



(evaluated at Q = Qo, A = 0) are negative or positive de- 
pending on the chosen direction. If we segregate absorb- 
ing (a) and conserving (c) nodes, the conductance matrix 
G may be put in block form. Two of them, the blocks 
G cc (real symmetric) and G ac will be relevant. This gives 
us the necessary tool to define the generalized minimal 
correlation coordinates. We consider the frequency de- 
pendent response by letting the evolution time extent to 
infinity, and introduce the time Fourier transform of the 
variables {Q c , A c , x°, where the vector {x c , X a } is a 
time-dependent source term introduced to produce cor- 
relation functions of the conserved and absorbed currents 

[see Eq. (O]. 

Following the steps of section HV Bl we again split the 
action into two parts, S — Sq + J dtV, where 



So 



, I dt[-A c Q c + A c G cc Q c + ( X c G cc + x a G ac )Q c 
dtvduj' 



2- 



[A c ( iW ' + G cc )Q c 
+ (x c G cc + x a Gac)Q c }6(oJ + uj'), (68) 



and where we have dropped the S in front of the variables 
for simplicity. As in Sec. II V Bl the generalized minimal 
correlation coordinates are defined by shifting and rescal- 
ing the A c variables in order to eliminate the x variables 
in Eq. (|68() . However, because \ is now a vector, the pro- 
portionality factor must be a frequency dependent ma- 
trix, 

A» - Dt(u;)[A» + GLx c H + <%,*»]. (69) 

Here D(uS) is the matrix network propagator, 

D(cj) = -(iuE + G cc )- 1 , (70) 

and E is the identity matrix. It is straightforward to ver- 
ify that after the shift, the functional J dtV becomes the 
generator of cumulants of minimal correlation currents, 
i.e. of the currents which are solutions of the Langevin 
equations: 



I c = 



)Q c a = -iu'Y^D a p{uj)I l 



(71a) 



0i 



I a a = ^Gc^DcrfiWifr+^icv, (71b) 

01a' 7 
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where I a p are the bare noise sources as defined in Eq. 
(J5J. We finally rescale X c ( w ) — * X c ( w )/(* w ) i n order to 
replace conserved currents with charges, I c — > Q c . 
The total action now acquires the following form 

S = {2m)- 1 J (LA c (-w)Q c (w) 

dtvW,D\A c + x c )+D ] Gl a x a ,X a ] , (72) 



where the simplified form of the A argument of V fol- 
lows after composing the various transformations. Fol- 
lowing the plan of the previous section, we replace the 
charge and counting variables {Q(uj) , A(oj)} by func- 
tional derivatives with respect to the charge and counting 



sources {J(w), K(u;)}, and take the V term outside of the 
functional integral. The functional integrals may now be 
performed to obtain 



W(3, K) = exp 



dojdoj' 
2tt 



j(t/)K(w)(5(w w; 



The perturbation V must now be expanded in a Tay- 
lor series with respect to all variables. The time depen- 
dence only appears through the variables themselves, so 
the expansion coefficients will be time independent, with 
the exception of the propagator D a p{u)) multiplying the 
counting variables. 



V 



gjiA hjn 



(Xl) 41 (Xr) 



I ! 



(((I^ ...(/; 



m i9 (Qi) ki ---(Q c P ) kp ))r 



v x (Qt) 31 my 

r*p ■ .11- Jn- 



(74) 



As in the one node case, the vertices 8q^{I%) vanish. 
We note again that the notation chosen for the expan- 
sion coefficients in Eq. I|74l) connects the formalism de- 
scribed here with the Langevin equation point of view. 
The minimal correlation cumulant ((. . .)) m may be cal- 
culated either by the expansion procedure described by 
Eqs. (|72I74|) . or by expressing the physical currents and 
charges in terms of the current source cumulants by solv- 
ing the Langevin equations for currents and charges, 
given by Eq. J7TJ ■ 

The nth order irreducible correlator 
((If (o;i) • • • Qn(u> n ))) ma y be expressed as a tree-level 
diagram with n external lines representing absorbed 
currents J„ and conserved charges Q c a . Every vertex 
is local in time, so if there are p legs at a vertex, 
each is assigned an independent frequency, while the 
time integral imposes overall frequency conservation, 
S(J2 i uJi). The cascade rules are generalized as follows: 

1. Every vertex represents the object 

$Qf(wi) ' ■■ S Qf(u l ){( I l+l(ui+l) ■ ••Qn( u} n)))m j 

which is multiplied by a (^-function conserving over- 
all frequency, <5(X)"=i ^i)- 

2. The minimal correlation cumulants 
((J; a +1 (w i+ i) ■ ■ ■Qn{oJ n ))) m may be evaluated by 
expressing them in terms of cumulants of sources 
{{I„a)) via the solutions (|7T)l of the Langevin 
equations, or by Eq. l|74|l if the Hamiltonian is 
known. 



3. The internal dashed arrow goes from Q c a (oj) to 
Sqc / u \ . It conserves the node index a and the fre- 

46 

quency u>. 

4. External lines for absorbed currents and conserved 
charges originate from I^{lo) or Q c a {io) of the ver- 
texes. They conserve the node index and the fre- 
quency. 

5. Sum over all internal node indices, and integrate 
over all internal frequencies to remove all but one 
of the frequency delta functions. 

6. The result has to be multiplied by the total number 
of inequivalent permutations. 

The cascade rules are easily extended to the field the- 
ory (see Sec. IHIfl . The functional analog to the inverse 
conductance matrix is the operator 

^(r-rO^^l^^^r-rOVW. (75) 

The diffusion propagator (iui + G -1 ) -1 can be used to 
solve the Langevin equations (|42|l for the density p(uj, r) 
and current I{u>) in order to evaluate minimal correlation 
cumulants. We would like to stress that these cumulants 
are limited to second order only, because in the diffusion 
limit the noise sources are Gaussian. The summation 
over node indices is replaced with an integration over the 
coordinate r. 
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V. APPLICATIONS 

The formalism presented above is intentionally ab- 
stract and general. This is to facilitate maximum appli- 
cability and not to tie it to a particular field. However, it 
is important to give concrete examples. For this reason, 
we give a detailed treatment of two problems. As a first 
problem, we consider the saddle-point equations of the 
ID field theories for D and F being arbitrary functions 
of the density p [see Eq. <|33[) ] . We apply the results of this 
analysis to the transport in a diffusive mesoscopic wire at 
zero temperature, rederive the FCS generating function 
of the transmitted charge obtained in Refs.Sandllil and 
give new results. We also prove the conjecture made in 
Ref. |52 that the current noise of the diffusive symmetric 
exclusion process at half-filling is Gaussian, i.e. all high- 
order cumulants of transmitted charge vanish. In the 
end of the Sec. IV Al we generalize our results to multi- 
dimensional diffusion models and prove the universality 
of their transport statistics. As a second problem, we 
address the statistics of charge fluctuations in a meso- 
scopic chaotic cavity. We explicitly find the probability 
distributions for different physical configurations. 



A. FCS for one-dimensional field theories. The 
mesoscopic diffusive wire. 

Before demonstrating our solution for the FCS of the 
mesoscopic diffusive wire specifically, we first consider 
the general ID field theory with the action l|33|) . In the 
stationary limit, p = A = the action can be written as 



1 



—Dp' A' + ^F(X') 2 



S = t J dz 
o 

The stationary saddle-point equations 



{FX' - Dp')' = 0, 2DX" + -f- (A') 2 = 0, 

dp 



(76) 



(77) 



can be partially integrated leading to the following two 
equations: 



Dp = ±y/l 2 - 2HF, 



X' = 2H/(X-Dp') 



(78a) 



(78b) 



The two integration constants X = —Dp' + FX' and 
H = -Dp'X' + (F/2)(X') 2 are the conserved (conditional) 
current and the Hamiltonian density, respectively. These 
conservation laws follow from the symmetries of our ID 
field theory [see Eqs. and lj39"|) and the surrounding 
discussion]. Thus we obtain the following result for the 
action lf7Bl . 



The equations l|78(l and l|79(l represent the formal solution 
of the FCS problem for ID diffusion models with D(p) 
and F(p) being arbitrary functions of p. The following 
procedure has to be done in order to obtain the cumulant 
generating function S(x) °f the transmitted charge: 



1. The differential equation l|78a|) has to be solved for 
p(z) with the boundary conditions p(z)\ z= o = pr 
and p(z)\ z= L = pr. The constant X should be ex- 
pressed through the constants pr, Pr, and TL. 



2. Next, p(z) is substituted into Eq. (|78b[) which is 
integrated to obtain A(z) with the boundary con- 
ditions Al = and Xr = x- 

3. Finally, using the solution for X(z) the constant Tt 
is expressed in terms of pr, Pr, X-, an d substituted 
into the action (|79|l . 

We note that by expressing H and x m terms of X, we 
may also formally obtain the logarithm of the current 
distribution, 



lnP(I)=S(I)-tl X (l), !->!, 



(80) 



as a result of the stationary phase approximation for 
the integral P(I) = J dxexp[S(x) — tlx] and because 
dH/dx = X/L. 

As an example of the ID field theory, we consider 
the FCS of the electron charge transmitted through the 
mesoscopic diffusive wire. When the chemical poten- 
tial difference Ap = pr — pr. > is applied to the 
wire, the electrons flow from the left lead to the right 
lead with the average current Iq — e~ 1 GAp, where G 
is the conductance of the wire. The elastic electron 
scattering causes non-equilibrium fluctuations of the cur- 
rent. At zero temperature, and for noninteracting elec- 
trons (the cold electron regime), the zero- frequency cur- 
rent noise power has been foundi^^^ to be equal to 
((I 2 )) = (l/3)e/o, i-e. the noise is suppressed compared 
to the Poissonian value. The suppression factor 1/3 was 
shown to be universalj22i£S i.e. it does not depend on the 
character of the disorder or on the shape of the wire. 
The FCS of the transmitted charge has been studied in 
Refs. l8l and lldl using quantum-mechanical methods, and 
recently in Ref. 22 using a classical method with the fol- 
lowing result for the generating function of cumulants of 
the dimensionless charge Q/e: 



S( X ) = (i/o/e) arcsinh 2 ^cxp(x) - 1 



(81) 



S = tLTL. 



(79) 



Here we will rederive this result using our classical 
method. 

On the classical level, the electrons in the diffusive wire 
are described by the distribution function f(z). Under 
transport conditions (and at zero temperature), this dis- 
tribution f(z) varies from Jr = 1 in the left lead to 
fit = in the right lead. Starting from the Langevin 
equation's as described in Sec. IHII or. alternatively, tak- 
ing the continuum limit for the series of mesoscopic 
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cavities^ we arrive at the action (|76(l with the for 
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1/2 



S={th/e) / dz[—f'X' + /(l — /)(A') S 



(82) 



-1/2 



where we have rescaled the coordinate z, p(z) has been 
replaced with the distribution f(z), and where D = 1, 
and P = 2/(1 — /) up to the overall constant Io/e. This 
form of F is quite general for fermionic systems. It orig- 
inates from the Pauli blocking factors, i.e. the transition 
probability is proportional to the probability that the ini- 
tial state is populated times to probability that the final 
state is empty^ Applying now the procedure described 
in the beginning of this section, we solve the saddle-point 
equations and find the fields / and A, 



/0,X) = g 



sinh(2az) 
sinh a 



X(z, x) = 2 arctanh [tanh(a/2) tanh(c 



a = arcsinh 



\J exp(x) - 1 



(83a) 



(83b) 



(83c) 



where TL = a 2 , so that according to the Eq. (|79|l we 
immediately obtain the result l|81|l . 

The logarithm of the current distribution ln[P(/)] can 
be now found from the equation l|8U|) . We obtain the 
following result: 



ln[P(7)] = -(i/ /e)[2acothaln(cosha) - a 2 



(84) 



where a has to be expressed in terms of X = I /Io by 
solving the equation 



a coth a = I/Iq. 



(85) 



The last equation has real positive solutions, < a < oo, 
for I > Iq, and pure imaginary solutions a = i[3 with 
< (3 < 7r/2, for I < Iq. The distribution P(Pj is 
strongly asymmetric around the average current I = Iq 
(see Fig. 0. It has the following asymptotics: InP = 
-(tl a /e)[l 2 - (2 In 2)1], for 1 = I/I » 1, i.e. P has a 
Gaussian tail, and InP = — (it 2 /4)(tl /e), for 1 = 0. 

We also plot the conditional electron occupation 
f(z,I), Eq. (|83a|l . for different values of the normal- 
ized current I/Iq. There are several interesting points 
to stress, (i) For large currents, I > Iq, the function / 
drops mostly at the ends of the wire, while for small cur- 
rents, I < Io, the drop of / is mostly concentrated in the 
center of the wire. This effect has a simple explanation. 
At the end points of the wire, z = ±1/2, the occupation 
f(z) is fixed independent of the particular value of the 
current I. On the other hand, its derivative takes the 
value /' = —X = —I/Iq at z = ±1/2, which can be eas- 
ily verified using Eqs. (|83a|l and l|85|l . As a result, f(z) 
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FIG. 7: The logarithm of the distribution of the current 
through a mesoscopic diffusive wire as a function of the ratio 
I/Io of the current to its average value Io. The distribution is 
strongly asymmetric, with the Gaussian tale at / 3> Io- Inset: 
The electron occupation / inside the wire as a function of the 
rescaled coordinate z, under the condition that the average 
current I = Io, no current 1 = 0, and large current I = 5Io 
has been measured. 



deviates from its linear behavior, f(z) = 1/2 — z, char- 
acteristic of the average value of current, 1=1. The ac- 
tual reason for this effect is that according to Eq. (|42b(l 
the total current 2 = —f' + is contains a contribution 
from the source of noise, v. The greatest contribution is 
concentrated at the center of the wire, where the noise 
power F = 2/(1 — /) has its maximum, while it van- 
ishes at the ends of the wire. Since the current X is 
conserved, /' has to be redistributed in such a way as 
to partially compensate the effect of the source v. (ii) 
Fluctuations of / are strongly suppressed at the ends of 
the wire, which is imposed by the boundary conditions, 
and at the center of the wire, as a result of the discrete 
symmetry, {z — > — z; f — > 1 — /}. (iii) Eq. (|8"5|l has ad- 
ditional solutions with (3 > ir/2. These solutions are not 
physical however, since / becomes negative or larger than 
1 leading to / < 0, which is impossible at T = 0. 

Returning to the saddle-point equations (j77(l . we note 
that if SF/5p = for a particular density po> then the 
fields p(z) = po, and X(z) = x z /L solve these equations. 
The fluctuations of the current become Gaussian with 
the noise power {{I 2 )) = F(po)/L. This generalizes and 
proves the conjecture made in Ref. |22| that the noise of 
the diffusive symmetric exclusion process is Gaussian at 
half-filling, /= 1/2. 

As a final remark we note that the whole class of multi- 
dimensional field theories, 



S = t [ dr[-VA£>Vp+(l/2)VAPVA], 



(86) 



with D = D(p) f,F = F(p) f, and f being an arbitrary 
constant symmetric tensor^ bear the same kind of the 
universality as the shot noise in diffusive conductors dis- 
cussed above (see Refs. ^ and l50|) . The reason is that 
the field theory with the action l|86|l can be mapped on 
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the ID theory with the action l|76(l by making use of the 
parameterization 

P (r)=pMr)}, A(r) = X[ip(r)], (87) 

where the function <p(r) satisfies the equation 

V • [TV<p(r)] = 0. (88) 

Using Eqs. 1)77(1 for p and A as functions of ip, it is 
straightforward to verify that the fields p(r) and A(r) 
given by l|57j) and satisfy the saddle-point equations 
for the action 1)86(1. One of the equations is the conser- 
vation of current: 

j = -DVp + FVX =lTVtp. (89) 

Since the ID Hamiltonian density is conserved, the action 
takes the following form 

s = tgn, (90) 

where the constant Q depends only on the geometry of 
the boundary d£l by Eq. JHHJ : 

9= drVcpfVcp = / ds-ipTVcp. (91) 
Jn Jon 

Consider now a two-terminal diffusive wire, so that the 
surface dfl consists of the left OVLl and right d^lji con- 
tact surfaces, and the open surface d£lo with no current 
through it. We choose the boundary conditions for tp to 
be 

<p{r)\ d n L = 0, <p(r)\ d n B = 1, ds ■ TV^(r)|an = 0, (92) 

so that p(r)\ d n L = Pl, p(r)\an R = Pr, K r )\an L = 0, 
AWIafijj = x, and ds ■ j(r)|an = °- Then TL becomes 
a function of p^, pn, and x, and the action 1(90(1 is the 
generator of the cumulants of the transmitted charge. If 
instead, TL and x are expressed in terms of X (as above for 
the ID theory), then one obtains the logarithm of the dis- 
tribution of the current, ln[P(J)] = S(I) — tlx(I), where 
I = QI, according to equations (JS5J), |J5TJ, and The 
constant Q may be interpreted as a "geometrical conduc- 
tance" of a wire. In particular, in the "ohmic" regime, 
i.e. when D is independent of p, we have Iq = X(x)\x=o = 
D(pl - Pr), and therefore Q = I /[D(p L - pr)}. In this 
case, the ratio 5/ Jo does not contain Q and becomes fully 
universal, proving also the universality of the result l(81|) 
as a special case. 

To summarize, we have proven the universality of the 
FCS of the transmitted charge for a two-terminal multi- 
dimensional generalized wire described by the action i(8t)|) 
with the noise tensor F(p)T ', being an arbitrary function 
of the charge density p, and with the constant diffusion 
tensor DT. The universality means that the FCS de- 
pends neither on the shape of the conductor, nor on its 
dimensionality. 53 The FCS of a mesoscopic wire given by 
Eq. 1(81(1 is a particular example of universal FCS. In the 
more general case, when D is a function of p. the FCS 
depends on the geometry through only one parameter Q, 
the geometrical conductance given by Eq. ((91|) . 




Q t /(NpAn) 

FIG. 8: The logarithm of the distribution of charge Q T in 
a symmetric cavity, Gl = Gr, averaged over measurement 
time r in the long time limit r ^> td and at zero tempera- 
ture. The results are presented for several transparencies F of 
the point contacts. It is clearly seen that the tails of the dis- 
tribution grow in the tunneling limit r<l. The distribution 
is symmetric, i.e. odd cumulants vanish. 



B. Charge fluctuations in a chaotic cavity. 

As another example of the applicability of the stochas- 
tic path integral approach, we now consider transport 
through a chaotic cavity. This problem is often investi- 
gated in mesoscopic physics because of its simplicity and 
conceptual clarity. A cavity consists of a large conduct- 
ing island of irregular shape that is connected to two 
metallic leads through quantum point contacts (see in- 
set of Fig. |SJ. The distinctive property of the chaotic 
cavity separating it from diffusive conductors is that the 
conductance is determined solely by the ballistic point 
contacts. The chaotic cavity itself may be either disor- 
dered or ballistic. Chaotic cavities can be described by 
a semiclassical theory if the point contacts have conduc- 
tances much larger than e 2 /h. The statistics of current 
flow through the cavity have been addressed using various 
methods. The zero-frequency noise power has been cal- 
culated using random matrix theory^ and the minimal 
correlation principle; 18 The higher order current cumu- 
lants have been obtained in Refs.0 andl20l The results 
are in complete agreement with random matrix theory. 

In this section we will address another type of statis- 
tics. In a typical experimental setup, the cavity is con- 
nected to the electrical circuit not only through the leads, 
but also through nearby metallic gates via the electro- 
static interaction. Observing potential fluctuations at 
these additional gates gives direct insight into the statis- 
tics of charge on the cavity. The noise power of the charge 
fluctuations in this system has been calculated in Ref.l55l 
The full statistics have been recently addressed using a 
random matrix theory^ Here, we rederive these results 
using the stochastic path integral, show new results on 
the temperature dependence of these statistics, and also 
investigate the instantaneous fluctuation statistics. 

In a semiclassical approach, both leads L, R and the 
cavity are described by electron distribution functions 
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/l>/r) and /. The Fermi functions in the leads f a = 
/f{E — fA a ) are characterized by their chemical potential 
\x a and their temperature T . The chaotic electron motion 
inside the cavity makes the cavity distribution function 
f(E,t) isotropic and position independent. Only its en- 
ergy dependence must be retained. From now on we set 
the electron charge to one, e = 1. Then the charge Q in 
the cavity is given in terms of the electron distribution 
function and density of states Np as Q = Np J dEf. 
The average value of charge is determined by the low- 
energy cut-off of the integral and is not relevant for the 
present discussion. The charge and electrostatic poten- 
tial of the cavity are related by a geometrical capacitance 
C g . In the following, we restrict ourselves to the case 
C g 3> e 2 Np which describes complete screening of the 
charge in the cavity. A more general discussion can be 
found in Ref. l56t To analyze the time evolution of the 
charge, we note that if the size of the cavity is smaller 
than the electron-electron and electron-phonon scatter- 
ing length, every electron entering the cavity at a certain 
energy leaves it at the same energy. The single electron 
energy is thus conserved and we can formulate a current 
conservation law separately for each energy interval dE, 



N F f(E,t) = J L (E,t) + J R (E,t), 



(93) 



where J a denote ingoing particle currents per energy in- 
terval dE in the left and right contacts. These currents 
are described by binomial processes with the cumulant 
generating function given by^ 

H a (f,iX a )dE = r-^dEln [1 +T/ Q (1 - /) (e lA - - l) 
+ r/(l-/ Q )(e- 4A °-l)], (94) 

where we have introduced the conductances of the point 
contacts G a , a = L,R, and their transparency T. 

The quantity of interest is the total number of electrons 
in the cavity averaged over the measurement time r, 



Qr = (N f /t) / dt / dEf{E,t). 



(95) 



We first consider the long time limit, r ^> 775, where Tp> — 
Nf/{Gl + Gr) is the average dwell time of an electron 
in the cavity. In this limit, the action is stationary with 
respect to the variables / and A, 

S = t fdE [H(f, iX) + i(N F /T) X f] , H = H L + H R , 



(96) 

where the external variable x generates the statistics of 
the desired quantity Q T . 

At zero temperature T = 0, the variables A and / are 
independent of the energy E, and the integration in Eq. 
(|96|l amounts to a multiplication by A/i = [Ll — fj, R . Eval- 
uating the Fourier transform of the characteristic func- 
tion Z(ix), 



(97) 



we express the full probability distribution P(Q T ) of 
charge on the cavity as an integral 

P{Q T ) = (27T)- 1 J dXexp[rAfiH(f,iX)] , 

/ = Qr/(N F Afi). (98) 

This integral will be calculated in the saddle-point ap- 
proximation. For the tunneling limit r <C 1 and for open 
point contacts T = 1 we obtain: 



lnP(Q T )=TGAfiK(f), 



K{f)\ 



r«i 



- fo) - Vfo(l - /) 



A'(/n ,., = /nln(£) +(l-/o)hi (J—j- 



(99a) 
(99b) 

, (99c) 



where G — Gl + Gr, and where we have introduced the 
average distribution function / = Gl/{Gl + G r ) in the 
cavity. We summarize that the results l|99(l have been 
obtained under the conditions T = 0, r ^> Td, and for 
r C 1 and r = 1. These results can be easily generalized 
to the case of a multi-terminal cavity. 

Although the general case of an arbitrary transparency 
r has been also solved analytically, the final expression 
for the charge distribution is too lengthy to be presented 
here. The Fig. [S] shows the distribution P(Q T ) at zero 
temperature for various transparencies T of the point 
contacts. The cavity is taken to be symmetric Gl = Gr. 
It is clearly seen that the tails of the distribution grow 
towards the tunneling limit. 

At finite temperature, further analytical progress can 
be made by considering the first few cumulants of the 
charge Q T . The integral l|9t)|) for the cumulant gener- 
ating function has to be evaluated at the saddle point. 
For x = the solution of the saddle-point equations 
dS/dX — and dS/df = are simply given by A = 
and / = f , where / = (G L f L + G R f R )/(G L + G R ) is 
the average electron distribution function in the cavity. 
From the diagrammatic technique discussed in Sec. II VI 
we derive analytical expressions for the first few cumu- 
lants. The second cumulant has been obtained in Ref.l55l 
As an example, we present here the result for the third 
cumulant for the case of open point contacts, T = 1: 



((Ql)) 



2t 3 d G l G r (G l - G R ) 
' r2 (G l + G R y 



F(An,T), (100a) 



F(Afj,,T) = A/i + 3 



A/i - k B T sinh(A^i/ k B T) 
cosh(A^/fc s T) - 1 



, (100b) 



Z(ix) = (27T)- 1 / dQdX exp(S) , 



where the function F(A[i, T) is always positive for A/j, > 
0. The first few cumulants are plotted in Fig. Has a 
function of the dimensionless bias Afi/k R T. Note that 
the fourth cumulant may change its sign as one goes from 
a symmetric cavity (/3 = 0) to an asymmetric cavity (/? = 
0.9). 

So far we have considered the time of measurement r 
longer than the dwell time tjj- Next we consider the op- 
posite limit r <C td (but still larger than tq — hj Afj.) and 
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FIG. 9: Cumulants of the charge inside a chaotic cavity, 
{{Qt))j n — 2,3,4 (in arbitrary units) as functions of the 
dimensionless potential difference A/i/fcsT. The parameter 
(3 = (Gl — Gr)/(Gl + Gr) characterizes the asymmetry of 
the cavity. 



study the instantaneous fluctuations of the charge Q in 
the cavity at zero temperature, T = 0. For this purpose 
we will use the stochastic path integral © for the propa- 
gator U(Qf, Qi, t) of the cavity charge. The distribution 
P(Q) of instantaneous fluctuations can be obtained by 
taking the t — > oo limit of the propagator U(Qf,Qi,t) 
and setting Qf — Q. We note that in the long time 
limit, (>to, the initial state Qi relaxes to the station- 
ary state Q, and as a result the saddle-point expression 
of the propagator U — exp(S' S p) factorizes according to 
S sp = So(Q) + Si(Qi) + Sf(Qf). Here the stationary con- 
tribution to the action is zero, So = 0, since there is no 
charge accumulation on a long time scale. We will show 
that the initial state contribution vanishes, Si = 0, so the 
system looses its memory about the initial state. Thus 
we obtain In P(Q) = S f {Q). 

We now focus on the case of a cavity with two tunneling 
contacts (r <C 1). Using the Hamiltonians in Eq. 1)940. 
and replacing the counting variable A — > i\, we write the 
action as 



S = GA/i / dt[r D \f + h s (X, /)] , (101a) 

K = (1 - /o)/(e A - 1) + /o(l - f)(e- X - 1) , (101b) 

where h s is the scaled Hamiltonian. The saddle point 
equations take the following form 

tdS = -(l-/o)/e A + /o(l-/)e- A , (102a) 
t d X = sinh(A) + (l-2/ )[cosh(A)-l]. (102b) 

The solution of the Eq. (|102b|> for A reads 

l + Af Q exp(t/T D ) 



X(t) = In 



1- A(l-/ )exp(t/T£>) 



(103) 



where A is the integration constant. 

To show that the initial contribution to the action Si 
is zero, we note that independent of the constant A, the 
absolute value of A is a growing function with the sta- 
tionary state given by A = at f = — oo. This means 



that starting from early times to — > — oo, the solutions 
areA(i) = and /(t)-/ = [f(t )-fo] exp[-(t-f )/r D ]. 
They describe the relaxation of the initial state f(to) to 
the stationary state / = /o- Substituting these solutions 
to Eqs. I|101|) we immediately find that Si = 0. 

After making this point we skip the rest of the details 
and present the final result for lnP(Q) = Sf(Q): 



In P(Q) | 



r<i = -TflGA/i 

+ (W) In 



i-/ 
l-/o 



(104) 



which can now be compared to the results (jHSJ- The 
cumulant generating function for the distribution l|104l) 
is given by S(x) = TnGA/i ln[l + f (e x - 1)]. Note 
that tdGA/i = A^pA/i is the total number of the semi- 
classical states in the cavity which participate in trans- 
port. Therefore the distribution l|104|l can be interpreted 
as being a result of uncorrelated binomial fluctuations of 
the Fermi occupations of each semi-classical state. We 
would like to mention that the same result can be ob- 
tained by solving the stationary master equation. 



VI. CONCLUSIONS 

We have put forth a stochastic path integral formu- 
lation of fluctuation statistics in networks. The math- 
ematical building blocks of the theory are 1) the prob- 
ability distributions of transport processes through the 
connectors, 2) a continuity equation linking the connec- 
tor currents to the charge accumulation in nodes (charge 
conservation) , and 3) a separation of time scales between 
nodal dynamics and connector fluctuations. The relevant 
action of the path integral is derived from these consider- 
ations and is related to the probability of (charge conserv- 
ing) paths in phase space. The dominant contribution to 
the statistics comes from the saddle point approxima- 
tion to the path integral, and the generating function for 
the interacting system is simply the action at the sad- 
dle point. Fluctuations are suppressed by the number of 
transporting elementary charges in the network. We have 
considered the continuum limit to obtain a field theory, 
and mapped it onto a Langevin equation with Gaussian 
noise. Cascade diagrammatic rules were found in agree- 
ment with Nagaev for the one node case, and extended 
to general current correlation functions in an arbitrary 
network. Applications to the current statistics of the dif- 
fusive wire and fluctuation statistics of the charge inside 
a mesoscopic cavity were also discussed. As the build- 
ing blocks of the theory are classical probability theory, 
the potential application of this formalism is very broad 
and applicable to any field where fluctuations are impor- 
tant, including mesoscopics, biology, economics, fluid and 
chemical dynamics. 

Note added in proof. — After this paper was submitted 
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for publication, the authors learned of previous related 
work by Bertini et alML Although they did not consider 
transport statistics, they did consider the probability to 
manifest a given macroscopic fluctuation of the particle 
density in diffusive lattice gas models and arrive at the 
action Eq. (J2HJ). However, the Gaussian nature of the 
local fluctuations was assumed a priori. We thank B. 
Derrida for bringing these papers to our attention. 
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